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Abstract 

We present a formal derivation of a simplified version of Compressible 
Primitive Equations (CPEs) for atmosphere modeling. They are obtained 
from 3-D compressible Navier-Stokes equations with an anisotropic vis- 
cous stress tensor where viscosity depends on the density. We then study 
the stability of the weak solutions of this model by using an intermedi- 
ate model, called model problem, which is more simple and practical, to 
achieve the main result. 

Keywords: Compressible primitive equations. Compressible viscous fluid, a 
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1 Introduction 

Among equations of geophysical fluid dynamics (see [8J), the equations govern- 
ing the motion of the atmosphere are the Primitive Equations (PEs). In the 
hierarchy of geophysical fluid dynamics models, they are situated between non 
hydrostatic models and shallow water models. They are obtained from the full 
3 dimensional set of Navier-Stokes equations for atmosphere modeling, 





(1) 


+ pdivV = 0, 


(2) 


D 1 D 


(3) 






(4) 


§11 = Q,^ 


RTp 


(5) 



where 

U is the three dimensional velocity vector with component u for horizontal 
velocity and v for the vertical one. The terms p, p, T, g stand for the density, 
the pressure, the temperature, the gravity vector {0,0, g). The diffusion term 
D is written as: 

D = mA.U + ud^yV (6) 

where A^^ stands for the derivatives of second order with respect to the horizontal 
variables x = {xi,X2), and p ^ v represents the anisotrope pair of viscosity. 
The diffusive term Qq represent the molecular diffusion where q is the amount 
of water in the air, and Qt is the heat diffusion standing for the solar heating 
(see for instance [17j for details of diffusive terms). The last term Cp is the 
specific heat of the air at constant pressure and R is the specific gas constant 
for the air. 

A scale analysis show that only the terms dyp and gp are dominant (see e.g. 
P,5j). This leads to replace the third equation of ([l} with the hydrostatic one to 
obtain the so-called Compressible Primitive Equations (CPEs) for atmosphere 
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modeling, 



^xP = D, 



dyp = -gp, 

—p + pdivU = 0, 
D ID 

§11 = Q,, 

RTp 

where x and y stand for the horizontal and vertical coordinate and 



(7) 



U • Va; + vdy 



Derivation "background" In this paper, we present the derivation of Com- 
pressible Primitive Equations (CPEs) close to Equations ([7]) (without taking in 
account complex phenomena such as the amount of water in the air and the 
solar heating) from the 3-D Navier-Stokes equations with an anisotropic viscous 
tensor. Emphasizing to the difference of sizes of the vertical and horizontal 
dimensions in the atmosphere (10 to 20 for height with respect to thousands of 
kilometers of length) , we derive the hydrostatic balance approximation for the 
vertical motion. We obtain simplified CPEs: 



dtp + div^ {p u) + dy (pv) = 0, 
dt {p u) -I- diva; {p-a.®u) + dy {p vvl) 
+dy {i^2dyu) , 
dyP{p) = ~9P 



^xP{p) = diva: {viD^{\i)) 



(8) 



where y stands for the vertical coordinate. The main difference between Model 
([7]) and Model (|8]) is the viscous term. Moreover, ii p = (? p with c = RT (for 
instance, as above), the density p is written ^(t, a;)e^^' ^ ^ where ^ , called again 
"density", is an unknown of the following system called model problem: 



dt 



9.e = 



(9) 



which is obtained from System (|S]) by the "simple" change of variables z = 1 — e ^ 
and w = e~^v where 

— = Si + u • Vj; + wdz 
dt 

and D stands for the following viscous terms 

D ^ div, {i^iD,{u)) + {iy2d,u) . (10) 
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As explained below, we cannot obtain a result directly on System ([H]) without 
using the intermediate Model (|9]), so we use the fact that System ([9]) is very 
close to System ([5]), and that the equation = is one of the key ingredient 
to achieve the stability of weak solutions of System ^ , to propagate the result 
to System ©. 



Mathematical "background" The mathematical study of PEs for atmo- 
sphere modeling were first studied by J.L. Lions, R. Temam and S. Wang (|12)) 
where they produced the mathematical formulation of System ([7]) in 2 and 3 
dimensions based on the works of J. Leray and they obtained the existence 
of weak solutions for all time (see also [17j where the result was proved by 
different means). For instance, in |17| . using the hydrostatic equation, they 
used the pressure p as vertical coordinate instead of the altitude y. Moreover, 
they wrote System (O in spherical coordinates {(j),9,p) to change compressible 
equations to incompressible ones to use the well-known results of incompress- 
ible theory. They distinguished the prognostic variables from the diagnostic 
variables, which are: (u, T, q) for the prognostic and (u, p, $) for the diag- 
nostic variables where $ is the geopotential gy{(l),9,p,t). Diagnostic variables 
{v — v{u),p = p{T),^ = ^{T)) can be written as a function of the prognostic 
variables through the div-free equation, p — RTp and by integrating the mass 
equation which is written in the new coordinates as follows: 

+ — = 0. 

P 

Then, the outline of their proof of the existence was: they wrote 



a weak formulation of the PEs (by defining appropriate space functions) 
of the form 

^ + At/ + B{U, U) + ElU) = I 
at 

where U — {u,T,q) with initial data U{0) — Uq and A,B,E are appro- 
priate functional, 

finite differences in time: J7", 

a priori estimates for C/", 

approximate functions: UAt(t) — on {{n — l)At,nAt) (following |16| ) 
a priori estimates for C/a*, 



and they proved the passage to the limit. 

Setting T and q constant, the main difference between Model (O and ((8)) 
comes from the viscous term © and ([TU)) . Starting from the Navier-Stokes 
equations with non constant density dependent viscosity and anisotropic viscous 
tensor, it is natural to get the viscous term ^TU^ . This term is also present in 
the viscous shallow water equations (see e.g. [fj). In the same spirit than |17| . 
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authors [H] showed a global weak existence for the 2-D version of model problem 
([8]) with p(p) = (? p (with the notation above) by a change of vertical coordinates 
(but not on p as done in [T7]) which led to prove that Model ([5]) and Model ^ 
are equivalent. Then, using an existence result provided by Gatapov et al |10| 
for Model they could conclude. Existence result [lOj for System © was 
obtained as follows: 

• by a useful change of variables in the Lagrangian coordinates, authors |10| 
showed that the density ^ is bounded from above and below . 

• by a "priori estimates and writing the system for the oscillatory part of the 
velocity u, authors [TO] obtained the existence result thanks to a Schauder 
fix point theorem. 

Unfortunately, this approach [TO] fails for the 3-D version ^ since the change 
of variables in Lagrangian coordinates does not provide enough information to 
bound the density ^. Moreover, to show a stability result for weak solutions for 
Model ([8]) with standard techniques also fails, since multiplying the conservation 
of the momentum equations of Model by (u, v) gives: 



d_ 

dt 



p\u\'^ + pin p - p + I dxdt + / i'i{p)\Dx{u)\'^ + nu2{p)\dyyu\dx+ / pgv dx 
n Jo, Jn 



where the sign of the integral / pgv dx is unknown (in the equation above, 

Jn 

\/ + V* u 

Dx{u) stands for ^ — ~)- It appears that, prima facie, there is a missing 

information on v to avoid the integral term / pgv dx introduced by the hydro- 

Jn 

static equation dyp = —gp. In fact, the study of the weak solutions stability 
cannot be performed directly on System (at least up to our knownledge): 
therefore, we have to study the intermediate Model ^ through the change of 
vertical coordinates. Indeed, we have just remarked that (as described above) 
that p is written as ^e~^l'^ where ^ does not depend on y. Thus, performing a 
change of variable in vertical coordinate in Model (I5|), following [9J, we showed 
that Model (|S|) could be written as Model As the hydrostatic equation in 
System ([9]) is dzt, = 0, an energy equality is easily obtained, and provided some 
a priori estimates. Nevertheless, those estimates are not strong enough to pass 
to the limit in the non linear terms; additional informations are required. On 
the other hand, we have to remark that the missing information for the verti- 
cal speed V for Model (or equivalently w for Model ([HI)) is fulfilled by the 
equation of the mass of System (|9]), which is also written as: 

dlzW = ^div^(^5yu) . 

Then, the fact that ^ — ^(t, x) combined with the equation above, allowed 
to obtain a mathematical entropy, the BD-entropy (initially introduced in [4], 
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where a simple proof was given in [TJ |B] or in [T] and the reference therein) . Let 
us also notice that, as for shallow water equations (see e.g. [2 [3l [5] to cite only 
a few), it is necessary to add a regularizing term (as capillarity of friction) to 
equations © (equivalently to Model ([S])) to conclude to the stability of weak 
solutions for Model in this present case, we add a quadratic friction source 
term which is written rpujul for System ([5]) or equivalently r^M|u| for System 
Indeed, the viscous term PU)) combined to the friction term brings some 
regularity on the density, which is required to pass to the limit in the non linear 
terms (e.g. for the term pu u, where typically a strong convergence of ^/pu 
is needed). Finally, by the reverse change of variables, the estimates, necessary 
to prove stability of weak solutions, were obtained for System ([5]) from those of 
System 

We note that, for the sake of simplicity, periodic conditions on the spatial 
horizontal domain fl^ are assumed, since it avoids an incoming boundary term 
(whose sign is unknown: see e.g. [B]), which appears when we seek a mathe- 
matical BD-entropy. Let us also precise that "good" boundary conditions on ft^ 
may be used (see |B]) instead of periodic ones to avoid this boundary term. 

We may also perform this analysis without the quadratic friction term by 
using the "new" multiplier introduced in |13| which provides another mathe- 
matical entropy: particularly to estimate bounds of pu^ in a better space than 



This paper is organized as follows In Section [21 starting from the 3-D 
compressible Navier-Stokes equations with an anisotropic viscous tensor, we 
formally derive the simplified Model ^ as described above. We present the 
main result in Section 12.21 In the third and last Section 13.21 we prove the main 
result. Firstly, we show that Model ([8]) can be rewritten as Model ^ which is 
more simpler. Then, taking advantages of the property of the density ^, adding 
a quadratic friction term (following [H [3]), we obtain a mathematical energy 
and entropy which provides enough estimates to pass to the limit in Model (|9]). 
Finally, following [S], the stability result for Model ([5]) is easily obtained. 

2 Formal derivation of the simplified atmosphere 
model 

We consider the Navier-Stokes model in a bounded three dimensional domain 
with periodic boundary conditions on fl^ and free conditions on the rest of the 
boundary. More exactly, we assume that motion of the medium occurs in a 
domain = {{x,y); x e 57^;, Q < y < h\ where ^l^ = is a torus. The full 
Navier-Stokes equation is written: 



L°°(0,T;L\f])). 



dtp + div(pw) = 0, 
df {pu) + div{pu (g) u) — drva{u) — pf = 0, 

P = P{p) 



(11) 
(12) 
(13) 
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where p is the density of the fluid and u = (u,z;)* stands for the fluid veloc- 
ity with u = (mi,M2) the horizontal component and v the vertical one. The 
pressure law is given by the equation of state: 

p{p) = (14) 

for some given constant c. The term / is the quadratic friction source term and 
the gravity strength is given as follows: 



f = -r^ul + ul{ui,U2,Q) -5k 

where r is a positive constant coefficient, g is the gravitational constant and 
k = (0,0, 1)* (where X* stands for the transpose of tensor X). The term 17(11) 
is a non symmetric stress with the following viscous tensor (see e.g. [TT | [T m [9]) 
E(p): 

y-iip) M2(p) 
M2(p) 

M3(p) M3(p) M3(p) 
The total stress tensor is written: 

ct(u) = -ph + 2S(p) : D{u) + A(p)div(u) h 

where the term S(/o) : D{u) is written: 



2^1 {p)D^ (u) /i2 (p) {dyU + Vxv) 
p.3{p) {dyU + \/ ^vf 2fi3{p)dyV 



(15) 



with It, the identity matrix. The term Dx{u) stands for the strain tensor, that 
is: j^, (u) ^ + where V.^^^-- 



dx2 



Remark 1 Let us remark that, if we play with the magnitude of viscosity fj,i, 
the matrix T,{p) will be useful to set a privileged flow direction. 

The last term A(p)div(u) is the classical normal stress tensor with A(p) the 
viscosity. The Navier-Stokes system is closed with the following boundary con- 
ditions on dH.: 

periodic conditions on dflx, 

v\y=o = v\y=h = 0, (16) 

dy-a\y=0 = 9yU^y=h = 0- 

We also assume that the distribution of the horizontal component of the velocity 
u and the density distribution are known at the initial time t — 0: 

u{0,x,y) = Uo{x,y), 

p(0,x,2/)=Co(x)e-^'/='^ ^ ' 

The fact that the initial condition for the density p has the form (|17p is justified 
at the end of Section [53] 

We assume that ^0 is a bounded positive function: 

< ^o(a;) < M < +00. (18) 
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2.1 Formal derivation of the simplified CPEs 

Taking advantages of the shallowness of the atmosphere, we assume that the 
characteristic scale for the altitude H is small with respect to the characteristic 
length L. So, the ratio of the vertical scale to the horizontal one is assumed 
small. In this context, we assume that the vertical movements and variations 
are very small compared to the horizontal ones, which justifies the following 
approximation. 

Let e be a "small" parameter such as: 

H V 

where V and U are respectively the characteristic scale of the vertical and 

horizontal velocity. We introduce the characteristic time T such as: T = — and 

the pressure unit P = 'pU^ where p is a characteristic density. Finally, we note 
the dimensionless quantities of time, space, fluid velocity, pressure, density and 
viscosities: 

~ t ^ X ^ y ^u_w 
t = — , X = y ~ — , u — — , V — — , 

P=^775, P=-, -^=T' /^j = = 1,2,3 
p X Mj" 

With these notations, the Froude number Fr, the Reynolds number associated to 
the viscosity pi (i=l,2,3), Rei, the Reynolds number associated to the viscosity 
A, Rex-, and the Mach number Ma are written respectively: 

F. = ^, i?e. = ^, iie. = ^, = (19) 

VgH Pi Ac 

Applying this scaling. System (fTT|)-(fT4l) is written: 

1 _ U V 

-d^p + — div£ (pu) + —dy {pv) = 0, 

pC/„ , , 'pU'^ , ~s pUV ^ , <?'P^ ~ 

—djXpu) + — ^divj {pu®u) ^ j^dy{pvu) + —^xP = 

^^^divjr {piD^{u)) + ^j^dyQI^dyu) + ^j-j^ dy (p^W ^v) + 

^ (a divj (2)) + AZ, (a dyv) , (20) 

pF„ , , pC/y ,. , , 'pV'^ ^ , OS c^fl„ _ 

—d^ipv) + — ^divj; {puv) + -j^dy [pv') + -jfdyp = 

-9PP+ ^JJf^^^x (M3 dyii) + ^^^divj (^ VjriT) + ^^^j^dy^pidyv) 

^TIl^^ pdiV£ (u)) + ^dy (XdyV^ . 

Using the definition of the dimensionless number (|19p , dropping ~ multiplying 
the mass equation of System (|20l) by T, the momentum equation for u of System 
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T T 
(I20p by — — , the momentum equation for v of System (120p by — — , we get the 
pU pV 

non-dimensional version of System (ITT|) - ((T4)) as foUows: 
dtp + div^r (p u) + dy (pv) = 0, 

dt {pu) + divj; {pu(»u) + dy (pvu) + —p^VxP = ^^divj; {piDx{u)) 

He I 

+ ^/12 ^^9yU+ Va;uj j + -^^"^x (Adiv:r(u) + XdyV) , 

dt (pv) + div, {puv)+ dy (pv^) + -^^dyp - -^-^P 

,2 p„ (Adiv^(u) + A9;,u) , 



(21) 



£2 

Next, if we assume the following asymptotic regime: 

^ = MP), ^ = sMp), ^ = 2, 3 and ^ = (22) 
and drop all terms of order 0(e), System (PT|) reduces to the following model: 
' dtp + divx (pu) + dy (pv) = 0, 



dt (pu) + div;i (pu® u) + dy {pvu) + -^\7 xp{p) = div^^ (i^iL'x(u)) 



+ dy {V2dy\l) + Pf, ^ ^ 

dyP{p) = 

called simplified CPEs. In the sequel, we simplify by setting Ma = F^. Then, 
the hydrostatic equation of System (|23l) with the pressure law (fH)) provides the 
density as 

p(i,x,2/) = e(i,a:)e-« (24) 

for some function ^ = C(^j2;) also called "density". Let us note that the density 
p is stratified: it means that for any altitude y, the density p has the profile of 
the function ^. Therefore, Equation ([M)) justifies the choice of the initial data 
(fT7|) for the density p at the time t = 0. In the sequel, we also assume that: 

y^ip) ^vp,i^ 1, 2, for v > 0. (25) 
2.2 The main result 

Assuming the viscosity under the form (P5)) and Ma — Fr, we define: 
Definition 1 A weak solution of System I12S\} on [0,T] x 57, with boundary Hid]) 
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and initial conditions HIT]) , is a collection of functions {p^u^w), if 
p e L°°(0, T; L'in)), VP e ^°"(0, T; H\n)), 

VVpeL2(0,T;(L2(r!))3) 

wzi/i p ^ and where {p, Vpu, ) satisfies: 

dtp + div^{y/p^u) + dy{^Uy/pv) = 0, 
p(0,a;) = po{x) 



(26) 



m </ie distribution sense, and the following equality holds for all smooth test 
function with compact support such as ip{T,x,y) = and ifQ = ipt=a- 

— / pudfLp dxdydt + / / {2vpDx{u) ~ pu® u) xfdxdydt 
Jo^Jn Jo ^ 

— / / pvvByip dxdydt — v / / pvByyip dxdydt + / / r p\u\iup dxdydt 
Jo^ Jn Jj) Jn Jo Jn 
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'iv{ip) dxdzdt + / / pv(pdxdzdt= / pottofo dxdy . 
Jn Jo Jn Jn 

(27) 

Now, we state the main result of this paper: 



Theorem 1 Let (/9„,m„,u„) be a sequence of weak solutions of System 
with boundary I116\} and initial conditions ilT]). satisfying entropy inequalities 
if 5' 71) and 1511) such as 



Pn>0, p'S ^ po mL\n), p^vJS ^ Po-uo ^ri L\n). (28) 
Then, up to a subsequence, 

• pn converges strongly in C°(0, T; L^^^{n)), 



/Ji^Ufi converges strongly in L'^(Q,T\L^^'^(VL) 



2\ 



• PnUn converges strongly in L^{Q,T;L^{VLf) for all T > 0, 

• {Pm \/PnUn, \/p7iVn) Converges to a weak solution of System 

• (pn,Un,Vn) Satisfy the entropy inequalities i f j'?)) and 151]) and converge to 
a weak solution of I12S\) - !!TS\) . 

The proof of the main resuh is divided into three parts: the first part consists 
in writing System (I23p . using (^, u,w — e^^v) as unknowns instead of (p, u, v). 
The obtained model is called model problem (see Section [XT]) . In the second 
part of the proof, we show the stability of weak solutions of the model problem 
(see Section [3.2.2113.2.^ . In the third and last part, by a simple criterion, the 
main result is proved (see Section 13.2.7^ . 
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3 Proof of the main result 



The first part of tlie proof of Theorem [T] consists in writing the simplified model 
l|23p in a more practical way, since the standard technique fails, as pointed out 
in Section [TJ 



3.1 A model problem; an intermediate model 

We first begin, by noticing that the structure of the density p defined as a 
tensorial product (see (1^^ ) suggests the following change of variables: 

z = l-e-y (29) 

where the vertical velocity, in the new coordinates, is: 

w{t, X, z) = e~^v{t, X, y). (30) 

Since the new vertical coordinate z is defined as — — z = e"**, multiplying by 

dy 

System ([^^ and using the viscosity profile and the change of variables (|29p 
provides the following model: 

9t^ + div, {in) + d,{^w) = Q, 

dt (C u) + divj, (^ u (g) u) + ^2 (^ u w) + ^ = vdlY, {iD,{n)) + vdU^n), 
d,X = 0. 

(31) 

which is the simplified CPEs with the unknowns 



{^{t, x),u{t, X, y),w{t, X, y)) instead of {p{t, x, y),u{t, x, y), v{t, x, y)) 

that we call model problem. In the new variables, the boundary conditions (|16p 
and the initial conditions ()17|) are written: 

periodic conditions on n^, 

W|2=o = w\^^h = 0, (32) 
dzU\^^o = ^^ui^^,, = 

and 

u(0,a;,y) = Uo(x,z), , , 

where n = T'^ x [0, h] with h=l- e"^ 



3.2 Mathematical study of the model problem 

This section is devoted to the study of stability of weak solutions of System (|3ip 
and equivalently for System as we will see in Section l3.2.7l In what follows, 
we can say that: 
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Definition 2 A weak solution of System on [0,T] x n, with boundary 
and initial conditions I133\) . is a collection of functions {^,u,w), if 



ith ^ and y^it, y^w) satisfies: 



e(0,x)=eo(x) 



(34) 



in the distribution sense, and the following equality holds for all smooth test 
function (/? with compact support such as (p{T, x, z) = and ipo — ipt=o: 



/ / S^vBtf dxdzdt + / / {2vS^Dx(u) — S^u® v) -.Vx^dxdzdt 

/ / S^wudz^ dxdzdt — V I I (^ud'^^^p dxdzdt + / / r£^\u\u(p dxdzdt 
jQ^Jn Jo Jn Jo Jn 

^divx{(p) dxdzdt = / ^oUoipodxdz. 
Jn 

(35) 

We then have the following result: 



Jfl 



Theorem 2 Let (^n, w„,w„) be a sequence of weak solutions of System iSl]) . 
with boundary Ii32\} and initial conditions ^33]). satisfying entropy inequalities 
and {HP such as 

6. > 0, a' ^ m LHn), ^ CoUQ i'(f^). (36) 

Then, up to a subsequence, 

• converges strongly in C^{0,T; L^^^{il)), 

• \/Cn'^n converges strongly in L^{0,T; L^^^{il)'^), 

• Cn'^n converges strongly in L^{0,T; L^{fl)^) for all T > 0, 

\/CnWn) converges to a weak solution of System 

• {^m UrnWn) satisfy the entropy inequalities I13T]) and 1151]) and converge to 
a weak solution of ll31]) - ![^} . 

The proof of Theorem [5] is divided into three steps: 

1. we first obtain suitable a priori bounds on u, w) (see Section [3. 2 .ip . 
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2. assuming the existence of sequences of weak solutions (^„,u„,w„), we 
show the compactness of sequences (^„, u„, Wn) in apropriate space func- 
tion (see Section inXMXni). 

3. using the obtained convergence, we show that we can pass to the hniit in 
all terms of System ([51]) : this finishes the proof of Theorem [2] (see Section 

MM- 

3.2.1 Energy and entropy estimates 

A part of a priori bounds on (^,u, w) are obtained by the physical energy in- 
equality which is obtained in a very classical way by multiplying the momentum 
equation by u, using the mass equation and integrating by parts. We obtain 
the following inequality: 

J [ (^y + (^lne-e + l))+ / a\DAn)\' + \d.u\')+r f ^\n\' ^ (37) 



which provides the uniform estimates: 

y/Cu is bounded in L°°(0, T; (L^ifl))^), (38) 

^/'u is bounded in ^^(o, T; {L^{n)f), (39) 

y^dzu is bounded in i^(0,T; {L^{n)f), (40) 

y^D^u) is bounded in 1^(0, T; {L^{n)f^), (41) 

Chi^ + 1 is bounded in L°°{0,T-, L^Q)). (42) 



As pointed out by several authors (see e.g. [31 [13]), the crucial point in the 
proof of the stability in these kind of models is to pass to the limit in the non 
linear term ig) u which requires the strong convergence of v/?u. So we need 
additional information, which may be for instance provided by the mathematical 
BD-entropy P]: 

we first take the gradient of the mass equation, then we multiply by and 
write the terms V^,^ as ^Vj; In^ to obtain: 

dt (2K V, In + div, (2i/^V, In e ® u) -t- a, (2i/^V, In ^w) 

+div^ (2KV^u) -f- i2iy^V^w) = 0. ^ ^ 

Next, we sum Equation (|43p with the momentum equation of System pip to 
get the equation: 

dt (e ^) + div, (^ ® eu) + d,{^w^)+ d, {2u^ V.w) , . 

+V,^ = 2i.div, (Mu)) - r^\u\u + v^d, (d^u) 

where tp = u + 2v\l ^ In^ and 2^2, (u) = V^u — V^u is the vorticity tensor. The 
mathematical BD-entropy is then obtained by multiplying the previous equation 
by '0 and integrating by parts: 
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{dt V) + div^ iij (E) Cu) + d,{^w V)) V' rfa;(iz = 37 / c^a^^^- 



(45) 

• Since 

/ 2iydivx{^Ax{u))\'xl'a^dxdz ^ 
Jn 

and periodic boundary conditions are assumed on fl^, we have: 

/ 2iydiv^ {£,A^{u)) u dxdz = / ^\A^{u)\'^ dxdz. (46) 
Jn Jn 

• Derivating the mass equation with respect to z provides the identity 

and also recalhng that ^ is only x-dependent, the integral becomes: 

/ dz{2iy^Vxw)Tpdxdz — / dz {2iy^V ^w) udxdz 
Jn Jn 

/ w d zdiv. {^u) dxdz (47) 
Jq 

1 



t\dzw\ dxdz. 

• The other terms are easily computed. We have: 

/ r^\u\uip dxdz = / r\u\^ dxdz + / 2iyr\u\uy ^S, dxdz , (48) 
Jo Jn Jn 

/ i^dz (dzu) ip dxdz = / i^S^ldzu]"^ dxdz . (49) 
Jn Jn 



and 



/ Wx^'ipdxdz = / W x£,u dxdz + 21^ / V^j^Va; In^ dxdz 
Jn Jn Jsi 

■ / ^+l)dxdz + 8iy / | V?l ^^^idz. 

Jn Jn 



7 

dt 



(50) 



Finally, gathering results (P5|) - ([50)l provides the following mathematical entropy 
equality: 

Jt /^(^^+^ln^-e + l)dxdz 

+ / i2i^^\dzw\'' + 2,^^\Axiu)\^ + i^^\dzu\^)dxdz (51) 
Jn 

+ / {r^\u\^ + 2iyr\u\u\/£, + 8i'\Wx\/£,f)dxdz ^0 
Jn 
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and estimates: 

V^/i is bounded in L°°(0, T; {L'^{n)f), (52) 

y/^d,w is bounded in L^{0, T; (£^(^2))), (53) 

yfyl:r(u) is bounded in 1^(0, T; (L^ {9.)f'^^) . (54) 

which finishes the first step of the proof of Theorem [2] 

Remark 2 Estimate is a straightforward consequence of estimates -y/fV" G 
L°°(0,T,L2(r2)2) and ^/lue L'^{0,T,L^{nf) since 

The second step of the proof of Theorem [5] will be divided into 4 parts. The 
first part consists to show the convergence of \f^i (see Section [3.2.2p . Then, we 
seek for bounds of \/^u„ and \f^Wn in an appropriate space (see Section [3.2.3p 
to be able to prove the convergence of ^„u„ (see Section [3. 2 .4^ . Thereafter, the 
convergence of \/^u„ is obtained (see Section l3.2.5p . 

3.2.2 Convergence of 

The first part of the proof of Theorem [2] consists to show the following conver- 
gence. 

Lemma 1 For every ^„ satisfying the mass equation of System {HP, we have: 
y/^ IS bounded m i°°(0, T, H^fl)), 

dt\/^ is bounded m L'^{0,T, H~^{n)). 

Then, up to a subsequence, the sequence ^„ converges almost everywhere and 
strongly in L^{0,T; L^{il)). Moreover, ^„ converges 

Proof of Lemma [1} 

The mass conservation equation gives 

IIv^(*)IIl2(o) = \\Q\\L^n)■ 



This equation and Estimate give the bound of in L°°{0,T,H^{n)). 
Using again the mass conservation equation, we write 

5t(V^) = -^v/?ndiv^(u„) - U„.Va;VC^- V^^z'^n 

^ r- ^ r- 

= 2 V CndlV:r(u„) - dlVa;(u„V^„) - V^„9zU;„. 

Then from Estimates (gT]), 

dty^ is bounded in L'^{0,T, H-'^{n)). 
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Next, Aubin's lemma provides compactness of \/Cn in C°(0, T, L^(r2), that is: 
y/C^ converges strongly to in C°(0, T, L^{Vt)). (55) 



We also have, by Sobolev embeddings, bomids of vsn in spaces L°°(0, T, U'{Vl)) 
for allpe [1,6]. 

Consequently, for p = 6, we get bounds of ^„ in £°°(0, T, Lr'{VL)) and we deduce 
that 

^„u„ = V^Vi^u„ is bounded in L°°(0, T, i3/2(f7)2)_ (55) 
It follows that dtS,n is bounded in L°°{Q,T,W''^'^'^{n)) since 

ftCn = -div(^„U„) - ^ndzWn 



and we have Estimate ([53| . 
To conclude, writing 

V,^„ = 2V6^V,V6^e i°°(0,r;i3/2(f^)2)^ 

we deduce bounds of ^„ in L°°(0,r; W^'^^'^{i} j). Then Aubin's lemma provides 
compactness of ^„ in the intermediate space 



compactness of in C"(0, T; ^^/^(fi)). 



3.2.3 Bounds of vCnUn and v^n^n 

To prove the convergence of the momentum equation, we have to control bounds 

of y^Un and y^Wn- 



Lemma 2 We have 

y^iLn bounded m L°°{0,T; {L^{n)f) 

and 

y/£,Wn bounded in L^{0,T; L^{il)). 



Proof of Lemma [2} We have already bounds of (see Estimates (pS)) ). 
There is left to show bounds of v^w^n ^^(0, T; L^(f2)). As ^„ = £,n{t,x) and 
Estimates (|53| holds, by the Poincare inequality, we have: 



h i-h 

2 J, ^ „ / IE) /■ rr„.. m2 



IvCnWnl dz^C \dz{y/ £,nWn)\ dz . 

Jo 



Hence, 



/ £,n\Wn\^ dxdz ^ C / ^„ jcJ^Wn |^ dxdz 

gives bounds of \f^Wn in L^(0, T; X^(fi)). 
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3.2.4 Convergence of ^„u„ 



As bounds of y^nUn and y' ^nWn are provided by Lemma [51 we are able to 
show the convergence of the momentum equation. 

Lemma 3 Let nin — £,nU,n be a sequence satisfying the momentum equation 
13 Then we have: 



^nthi ->-m in L'^{0,T; {LP{n)f) strong ,V1 p < 3/2 

and 

£,nUn m a.e. {t,x,y) € (0,T) x fl. 

Proof of Lemma [3) 

Writing Va;(CnU„) as: 

provides 

V,(e„u„) bounded in L^O.T; (L^n))^''^). (57) 

Next, we have 

d.i^nUn) - v/e^vCa,(u„) is bounded L2(0,r; {L^/^fl))^). (58) 
Then, from bounds (|57|) and (|58)) . we deduce: 

^„u„ is bounded L^{0, T; W^'\nf). (59) 
On the other hand, we have: 

dt{£.n U„) = -div^; ($„ U„ (g) U„) - (^„ U„ Wn) - Va;Cn 
+ i^diVa; (Cn£'x(u„)) + (C«9zU„) . 

As 

f„u„ «) u„ = v/?u„ «) v^u„, (60) 

we deduce bounds of 

e„u„(»u„ in L==°(0,T;(Li(f^))2x2). 

Particularly, we have 

div(^„u„ «) u„) bounded in L°°(0,T; (l^-2.4/3(-f^))2^^ 

Similarly, as ^nUnUi^ = ^ (^^(^))^: ^'l^'^ have 

5.(e„u„«;„) bounded in L°°(0,r; (W^-2,4/3(^))2)_ 

Moreover, as 

&.V^d,u„,e L^{0,T; {L^/\n)f) and 
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we get bounds of 

We also have bounds of V:^^„ G L°°(0,T, (T/F"^'^/^(17))2). 
Using T^"i'3/2(^) ^ W-^'^/'^iVt), we obtain 

9t(^nU„) bounded in L^{Q,T;W-^''^'^{Vlf). (61) 

Using bounds ([5^ . (pT|) with Aubin's lemma provides compactness of 

e„u„ e L\G, T- iLP{n)f), yp e [l, 3/2[. (6^ 

3.2.5 Convergence of \/^u„ and ^„w„ 

Let us note that, up to Section |3.2.4[ we can always define u — m/^ on the set 
> 0}; but we do not know, a priori, if m equals zero on the vacuum set. To 
this end, we need to prove the following lemma: 

Lemma 4 

1. The sequence y/^Un satisfies 

• VCn'^n converges strongly in L^(0,T; L^{n)) to —=. 

• We have m — almost everywhere on the set — 0} and there exists 
a function u such that n2 = S^u and 

£.nUn £,u strongly in £^(0, T; LP{n)'^) for all p € [1, 3/2[, (63) 



V^iUn ^/lu strongly in L^(0, T; L^{nf). (64) 
2. The sequence \/^Wn converges weakly in L^{0,T; L^{il)) to \/£,w. 
Proof of Lemma |4j 

We refer to [14J for details of the first part of the proof. The second part of the 
theorem is done by weak compactness. As is bounded in L^{0, T; L^(i7)), 

there exists, up to a subsequence, \/^w„ which converges weakly some limit I 
in L'^{0,T;L'^{n)). Next, we define w 

[ a.e. if ^ = 

where the limit I is written: I = y^— ;= = \/S,w. 

V? ■ 

This finishes the second step of the proof of Theorem |5] 

In the third and last step (see Section I3.2.6P , using the convergence of Sec- 
tions 13.2.214372751 we show that we can pass to the limit for all terms of System 
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3.2.6 Convergence step 

We are now ready to prove that we can pass to the hmit in all terms of System 
(PT|) in the sense of Theorem[5] To this end, let (^„, u„, w„) be a weak solution 
of System ([5T|) satisfying Lemma [T] to |3] and let e C^([0,r] x fi) be a smooth 
function with compact support such as (/'(T, x, z) = 0. Then, writing each term 
of the weak formulation of System ([5T|) . we have: 



/ / dt{S,nUn)4>dxdzdt — / £^nUndt(j) dxdzdt 

Jo Jn Jq. Jn 

£,QUQ(f>{0,x,z)dxdz. 



(65) 



n 



Using convergences (|36l) and Lemma [31 we get 

— / / ^„Undt(f> dxdzdt ~ / ^qUq 0(0, x, z) dxdz —> 
^0 Jn Jn 



/ ^udt4> dxdzdt ~ / £,oUo(j){0,x,y) dxdz. 
Jn Jn 



/ / div X {£,n'^n 'S) Un) ■ (f> dxdz dt = — / / ^„u„ (g) u„ : Va;(/) dxdzdt. 
Jo Jn Jo Jn 

From Equality and Lemma 01 we have: 

/ $nU„ (X) u„ : Vaj^dxdzdt — > — / / (8) u : Vaj^dxdzdt. 
Jn Jo Jn 



/ / dz{£,nUnWn) ' (f) dxdzdt = - / / C„u„w„ • d^cfidxdzdt. 
Jo Jn Jo Jn 

As ^„u„w„ = \/l^\in\/^iWn, by Lemma m we get: 

/ ^nU„?i;„ • dz4> dxdzdt — / ^uw • dz4> dxdzdt. 

/ / ^ x^7i ■ (f> dxdzdt — — / / ^ndiva: ((/)) (ia;(iz(ii. 
^0 Jn Jo Jn 

Then, Lemma [T] provides: 

^ '■ . r f . 

^ndivx{(f>) dxdzdt ^ — / / ^divx{(f>) dxdzdt 



"'n "'0 Jn 
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/ / divxiS.nDxi'Hn)) ■ <l)dxdzdt ^ - / / £,nDx{u.n) : Vcjidxdzdt. 
Jo Jn Jo Jn 

Since Dx{un) — -(Va:U„ + V^u„), expanding the term in the last integral 



gives: 



/ / ^nDxiun) ■■ W xcj) dxdzdt 
Jo Jfl 

O / / i^n^n ■ ^x<t> + ^ x4>^ xiVin) ' \/TnVin dxdzdt 

- / / (^„u„ • diva;(V^(/)) + V^y^n • ■ \JTn^n) dxdzdt. 
^ Jo Jn 



From Estimates ([52]) . the sequence \7x\/S,n weakly converges, and using 
Lemma [U Lemma [3] and 21 we obtain: 



T; [ [ (CnUn • '^x<l> + ^x4'^x{\/Jn) ' a/^U" dxdzdt 

^ Jo Jn 

\ I I (C«u„ • div^(V^0) + V^v^ • • v^Un) dxdzdt 
^ Jo Jn 

\\ f {S.U- Ax(j) + Vx(j)Vx{y^)- y/^u dxdzdt 
^ Jo Jn 

+ 11 I {£.u ■ div^(V^0) + V^y? • "^xcf^ ■ V^u) dxdzdt. 
^ Jo Jn 



Hence 

i-T 



\ \ £,nDxiun) :Vx4> dxdzdt — / / S,Dx{u) : V ^^j) dxdzdt. 
Jo Jn Jo Jn 

/ / dl^{^nUn) ■ (j) dxdzdt ^ / / ^„Un ■ d^^{4>) dxdzdt. 
Jo Jn Jo Jn 



lo Jn Jo Jn 

Using Lemma [3] provides the following convergence: 



T 

2 



Jn 



^„u„ • d^^{4>) dxdzdt 



I / S^vl ■ dl^{(l)) dxdzdt 
Jn 



i i rS,n\vin\'Vi.n ■ 4> dxdzdt ^ I I rS\\i\\i ■ (f) dxdzdt 
Jo Jn Jo Jn 

with Lemma |4j which finishes the proof of Theorem [2] 
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3.2.7 End of the proof of Theorem [T] 

In order to conlude, let ((^„, u„, w„) be a weak solution of System ([5T|) . then 
all estimates 13. 2. 211X2. 61 hold if we replace ^„ by p„ and w„ by w„ (see [S]), 

since p(t,x,y) = ^(i:,x)e~^ and w{t,x,z) ~ v(t,xy)e^'^ where —z — . This 

ay 

proves Theorem [TJ 
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